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ABSTRACT 



We study the effect of tlie dynamicaf friction induced by tlie presence of 
substructure on tlie statistics of the collapse of density peaks. Applying the 
results of a former paper (Antonuccio-Delogu and Colafrancesco, 1994) we 
show that within high density environments, like rich clusters of galaxies, the 
collapse of the low v peaks is strongly delayed until very late epochs. A bias of 
dynamical nature thus naturally arises because high density peaks preferentially 
collapse to form halos within which visible objects eventualy will form. We 
then derive a selection function for these collapsing structures: with this in 
hands, we can calculate the values of the bias coefficient on cluster scales for 
any hierarchical clustering scenario. For a standard CDM model we show here 
that the dynamical bias that we derive can account for a substantial part of the 
total bias required by observations on cluster scales. 



Subject headings: cosmology: theory — galaxies: clustering — galaxies: 
formation 
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1. Introduction 

Several observational evidences indicate that the distribution of groups and clusters of 
galaxies has to be biased with respect to the underlying distribution of cosmic material on 
large scales. A widespread theoretical motivation for this is to assume that the formation 
of actual structures proceeds around the high peaks of the initial density field (Kaiser 1984, 
Peacock & Heavens 1985, Bardeen et al. 1986 - hereafter BBKS). 

Within this frame (and assuming a Gaussian statistics) the peak correlation ^pk{r) ~ h'^ip{r) 
is directly related to the structure of the underlying density field: hence, the power and 
mass spectra are proportional through the biasing factor h\ {6'^{k)) = 6^(A^^gg(A;)). This 
Lagrangian biasing is spatially uniform; it intervenes also in the normalization of the overall 
amplitude of the linear perturbation spectrum through the function J3(r) = / dr^[r)r^ . In 
terms of this Lagrangian biasing, the 2-point correlation function for rich clusters indicates 
6 ~ 3 (see, e.g., BBKS), a result also confirmed by the APM survey (Maddox et al. 1990). 
Moreover, the observed mass distribution of Abell clusters (Bahcall & Cen 1992) requires 
a strong suppression of galaxy formation. Within the spherical model for the evolution 
of density perturbations (see Peebles 1980) this is consistent with a value b 1.7 1.9 
(see also Antonuccio-Delogu and Colafrancesco, 1994). An extensive analysis of the X-ray 
luminosity functions and of the X-ray counts for galaxy clusters (see Colafrancesco and 
Vittorio 1994, hereafter CV) indicates that the local and distant distributions of X-ray 
clusters depend on the combination b ■ 6c ^ where 6c is the linearly extrapolated threshold 
for collapse. Values 6 • (5c ~ 3 obtain for a wide range of cosmological models. This result 
is also confirmed by recent hydrodynamic simulations (Cen & Ostriker 1993, Bryan et al. 
1994) showing that a COBE normalized (unbiased) standard CDM scenario predicts a 
space density of X-ray clusters in excess by a factor ~ 5 with respect to the observed value. 

From the result of CV, it appears that the cluster abundance depends on the effects 
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of non-linear dynamics during cluster formation (for which 6c is an indicative parameter) 
and/or from the uncertainty in the normalization of the linear power spectrum, as contained 
in the value of the bias parameter b. Both these issues are emanations of a biasing process 
already effective. 

The physical origin of such a biasing is not yet clear even though several mechanisms 
have been proposed (see Rees f985, Dekel & Rees f987, Dekel & Silk f986, Carlberg 
f99f , Cen & Ostriker f992, Bower et al. f993, Silk & Wyse f993). But, if efficient, these 
mechanisms will generate a physical selection of those peaks in the initial density field that 
eventually will give rise to the observed cosmic structures. 

In this paper we will derive (in §2) a selection function for the peaks that give rise to 
the protostructures: it results from a proper treatment of the dynamical effects of small 
scale substructure during the early stages of the collapse. In §3 we will derive values of the 
biasing parameter on the relevant filtering scales as produced by the previous non-linear 
mechanisms and we will discuss their cosmological relevance in the conclusions (§4). 



Following BBKS we define a selection function t[p) as the probability that a peak of 
the filtered density field will eventually form an object: 



where Uobj^ Npk{i') are the objects' and peaks' densities, respectively (we restrict ourselves 
to overdense regions, so the lower limit of integration in eq. I is put equal to 0). The 
specification of the selection function is crucial to relate the properties of the "objects" 
forming out of the density peaks. BBKS adopted an empirical selection function of the 
form: t[p) = [v jvtY l[l -\- [v jvtY], which depends on two parameters, namely the threshold 
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Selection function 
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and a shape exponent q. The results they obtain concerning the statistical properties 
of clusters and galaxies depend then on these parameters, because the selection function 
describes the statistical relationship among peaks and objects. 

We will now give a more physically motivated prescription to obtain a function t[p) selecting 
precisely the "objects" we defined as galaxy clusters. Let us consider the average density 
profile around a peak of a given central height p (see Figure 1) in the density field smoothed 
on an arbitrary scale Rj. The outermost shell defining our objects will be at a distance 
defined in such a way that the total mass contained within the average density profile of 
the peak is fixed to an arbitrary value Mt- This distance depends also on p and on and 
thus in the following we will compute the bias coefficients for different values of Mt and Rj. 
The average overdensity at the distance can be either larger or smaller than the critical 
overdensity for collapse. In the latter case, if the density profile of our peak coincides with 
the average one, the shell will not collapse. However, the previous argument holds for the 
average profile: at any distance there is a finite probability that the peak overdensity will 
be different from the average one, and for a Gaussian density field this probability is still a 
Gaussian: 

p(MS>M)=^exp(-^^|Ki£) (2) 



where {S){r) is the average overdensity and the dispersion cr§(r) is given, e.g, by eq. 3.8 
of Lilje & Lahav 1 99 1. Given a critical overdensity threshold for collapse (5c, for any given 
density peak of specified height p we can define a probability that the outermost shell will 
collapse as: 



t{p)= / dpp[6,{6){rMr)\ (3) 

The integrand in eq.(3) is evaluated at a fixed distance vm^- We define t[p) in eq.(3) as our 
selection function. Its dependence on p and Rf is implicit through the dependence from 
the average density profile 6[r). For p ^ oo the average density raises much over (5c, so one 
has 6c <^ '5(rMT), and t[p) — )• I. In the opposite case, when ;/ — )• 0, the variance (t| in eq. 2 
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is always very large over the integration domain, so t[p) — )• 0. In Figure 2 we plot t[p) for 
different values of the lower integration limit in eq. 3. 

As we have shown in a former paper (Antonuccio-Delogu and Colafrancesco, 1994, hereafter 
AC94), the presence of small-scale substructure modifies the value of the critical threshold 
for collapse 6c from the value Sco ~ 1.68 holding for the ideal case of a collapse without 
substructure. The dynamical friction induced on objects contained within infalling shells 
by the small scale substructure induces a drag force which act to contrast the gravitational 
force. As a consequence, the value of 6c is larger w.r.t. the ideal case (AC94): 

Ao 



Sc = S. 



cO 



1 + 



I - ij.{6c 

where the dynamical friction coefficient is given by (see AC94, eq. 23): 



(4) 



Ao = 4.44^M^ . log 

{V}av 



I.I2- 



1/2 



TT 



SttGpc 

and the function fi[6) is given by: 



(5) 



In eq. 5 the quantities in bracket are averages over the population of "subpeaks" generating 
the small scale substructure. The quantity c[6) was defined in AC94. 

Eq. 4 is an implicit equation for 6c which can be solved once Aq has been specified. This 
latter quantity depends on the mass spectrum and temporal evolution of the small scale 
substructure. These two important inputs are today poorly known both from theoretical 
and from observational grounds. We then decided to adopt the same strategy as in AC94, 
namely to include only the contributions from the peaks selected up to a maximum mass 
of 7.34 X 10® Af©, which is more than three orders of magnitude below the filtering mass 
of the considered objects (hereafter, we will refer to this population of density peaks as 
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the small-scale substructure , SSS). The values of Aq so obtained is then a lower limit, also 
because we do not take into account in eq. 5 any term depending on the correlation function 
(Antonuccio-Delogu, 1992;Antonuccio-Delogu and Atrio-Barandela, 1992). However, 
we take into account the fact that the peaks' number density will be enhanced within 
protoclusters w.r.t the background. Following BBKS the enhancement factor is given by: 



/e 



exp 



(7) 



^ ohj ^ \ ^obj J 

where a subscript "c/" denotes quantities computed within the protocluster. In the following 
of this Letter we will show results obtained for Pci = 2 and = 5h~^ Mpc, typical of a rich 
Abell cluster. 



3. Bias coefficients 



The bias coefficient of a given class of objects can be given by (see BBKS): 



{Rf) = ^ + 1 

where the average {p) is given by (Bardeen et al. 1986, eq. 6.45): 



;t{p)N.pk{iy) 



(9) 



1 -72 

and 7 and 6 denote the usual spectral quantities (see BBKS and AC94). Once we have 
specified a given spectrum we have all the tools to evaluate the integral in eq. 9 and then 
to calculate the bias parameter bss- 

For the sake of brevity, we will now restrict ourselves to the standard CDM model 
leaving to a more general paper the applications to other scenarios of structure formation 
(Antonuccio-Delogu, Colafrancesco and Vittorio 1994). At present, there is a rather large 
uncertainty about the value of the normalization constant of the standard CDM spectrum 
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with values ranging from 63 ~ 1 (coming from the analysis of COBE data, Scaramella and 
Vittorio 1993) up to bg ~ 1.8. Here we present our results for four different values of bg^ 
namely bg = 1.2, 1.4, 1.6, 1.8, and for three values of the Gaussian filtering radius of 
the objects: Rj = 200,356 and 500 kpc. For the first two values of Rj we choose 
in such a way that the total mass enclosed by the outermost shell defining the object is 

2 X lO^^Af©, a value typical of a bright spiral galaxy (see Taable 1). So from a comparison 
of the results for these two cases we can have insights on how the bias induced by small 
scale substructure is affected by the compactness of the structure. The results are shown 
in Table 1 and Figure 3. In the fifth column of Table 1 we report the values of bss^ the 
dynamical friction-induced bias coefficient, computed from eq. 8 and from the proper values 
of 6c (see eq 4). For larger Rj (at fixed Mt) the value of bss increases: this is a consequence 
of the behaviour of the selection function t[p) in these two cases. In fact, t[p) acts as a 
filter in eq. 8, and from Figure 2 we see that the approximate value of p at which t[p) ~ 1 
increases with increasing Rj. The physical reason for this is that, leaving constant the value 
of Afx, a peak selected with a larger filter radius is shallower than a peak selected with a 
smaller radius, and the gravitational force on the outer shells will be smaller. 
The error bars in Figure 3 were computed by taking into account the theoretical uncertainties 
in the calculation of Aq. These originate from the fact the quantities {m)av^ {n)av depend 
on very uncertain functions like the typical mass associated with a peak of a given height, 
mpeak{v). We have adopted the formula given by Peacock and Heavens (1985), but in order 
to compute the error bars we have also allowed the value of mppak{v) to be estimated by a 
top-hat filtering mass. Analogously, we have taken into account the iblcr on the number 
density of the peaks (BBKS). 

We also notice that for smaller peaks (i?/ = 200/i"^ Kpc, Mj = 1.4 X lO^^M®) the 
bias coefficient is also smaller, consistently with the interpretation we have just given. In 
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general, the values of bss are smaller than the global bias coefficient bg^ as it should be. 
However, they can account for a large fraction of the total biasing level, particularly at low 
values of bg which seem to be favoured by COBE. 

4. Conclusions 

In this Letter we have looked at one of the consequences of the delay of the collapse 
of protostructures induced by the presence of small-scale, fine-grained substructure. 
This collapse delay requires that the critical overdensity for the collapse of a spherical 
perturbation has to be larger than the canonical value 6c > Sco ~ 1-68 (see also Antonuccio- 
Delogu and Colafrancesco, 1994). 

As a consequence of this result, within overdense regions on the cluster scales, the 
progenitors of the structures which collapse and virialise by the present time are among the 
peaks which are selected according to the selection function t[p) introduced in section 3. 
This statistical selection function takes into account the dynamical biasing of realistic 
protostructures out of the structure of the density field at early times. Thus it is not a 
pure Heaviside function as in the empirical case treated by BBKS: its shape depends on the 
assumed overdensity threshold for collapse 6c and on the probability distribution (see eq.2) 
for having a collapsing shell around a density peak. Moreover this function selects peaks of 
larger height for increasing filtering scale Rj. 

Based on our selection function, we have also derived detailed values of the dynamical 
bias, bss, on the cluster scales: we found that bss ^ ^8- Depending on the values of 6c and 

we found that bss/bg ~ 0.5 0.8 for values bg > 1.2 as indicated by the bulk of the 
observations on cluster scales. This means that a sensitive fraction of the bias coefficient bg 
required to match the observations of the cluster abundance can be accounted for by the 
non-linear, dynamical effects of substructure during the formation of galaxy clusters. 
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Here we have also shown that the non-linear effects during cluster collapse producing 
large values of 6c traslate consistently in high values of the biasing parameter bss > 1- This 
justifies theoretically the conclusions of CV, although those were based only on a linear 
approach. 

These results are obtained specifically for a standard CDM scenario (Oo = 1, /i = 0.5): 
however, they hold in principle for any hierarchical clustering scenario. Their detailed 
importance in different models of structure formation will be discussed in another paper. 
In fact, a detailed calibration of the importance of this and other non-linear mechanisms 
during cluster formation and evolution is needed to assess definitely the origin of the 
cosmological biasing in the theory of structure formation. 

Acknowledgments. The investigation of V.A.-D. was partly supported by Danmarks 
Grundforskningsfond through its support for an establishment of the Theoretical 
Astrophysics Center (TAG). 
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Table 1: Bias coefficients. 

"In h-'^ Mpc 

^In units of W^Ma, 
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Figure Captions. 



Figure l.-The average density profile of a peak [continuous curve) together with the zblcr^ 
profiles [dashed curves). The asterisk marks the lower integration limit of the integral in 
eq. 3 (see section 2.). 

Figure 2. -The selection function t[p) for Rj = 356 Kpch~^ [6c = Sco [dots), 6c = 3.51 
(long-dashed)) and for Rj = 500 Kpch~^ [6c = 6co [short-dash), 6c = 3.51 continuous)). 

Figure 3. -Bias coefficients. We plot the values of bss for Rj = 500h~^ kpc. Filled circles: 

= 6df, open circles: 6c = 6co. . The dotted line marks the maximum value of bss (i-e. 
bss = bg). The error bars have been computed by summing the squares of the uncertainties 
in bss due to independent uncertainties on (rn)av, {n)av and on the central correlation 
function i^o (Antonuccio-Delogu, 1992). 



